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We show that using the modified form of the Dirac Hamiltonian as suggested by Bekenstein does 
not affect the analysis in QL |3 U 9 OI QSO data pertaining to a measurement of a variation. 
We obtain the present time limit on Bekenstein's parameter, tan 2 \ = (0-2 ± 0.7) x I0 -6 , from the 
measurement of the hydrogen 2p fine structure using value of a obtained from different experiments. 
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I. INTRODUCTION 

The possible variation of the fine structure constant, 
a, is currently a very popular research topic. Webb et 
al. have found evidence of a variation by analyzing 
absorption lines in QSO spectra, while other groups 
have used the same method [4| but found no evidence of 
a variation. Recently Bekenstein p| has questioned the 
validity of the analysis used by these groups. Bekenstein 
shows that within the framework of dynamical a vari- 
ability the form of the Dirac Hamiltonian relevant for 
an electron in an atom departs from the standard form. 
Unfortunately no self consistent quantum electrodynamic 
theory was derived. Instead the Dirac Hamiltonian, H, 
was presented for an electron bound by a Coulomb field: 



H = H + SH 
H = (—ihcot ■ V + mc 2 f3 - 
SH = (I - (3) tan 2 X • V c 



e$I) 



(1) 

(2) 
(3) 



where Vc = —Ze 2 /r. The last term is related to an effec- 
tive correction to the Coulomb field due to the dynamical 
nature of a, and tan 2 \ is a small parameter. In other 
words, a dynamically varying fine structure constant can 
be accounted for as a perturbation of the Dirac Hamil- 
tonian. This perturbative term, SH, vanishes in a non- 
relativistic approximation but produces some relativistic 
corrections which can be studied both in astrophysical 
spectra and laboratory conditions. 

In the following sections we will show how this pertur- 
bation shifts atomic energy levels. In particular we pay 
attention to heavy atoms which provide us with astro- 
physical data and are the most sensitive to a possible a 
variation. We also consider atomic hydrogen, since it is 
the best understood atomic system for laboratory exper- 
iments. 



II. MULTIELECTRON ATOMS 

Multielectron atoms are of interest to us since they 
are the most sensitive to a varying a. We performed a 



calculation to show how the modified form of the Dirac 
Hamiltonian affects the energy of an external electron in 
a heavy atom. We averaged (5H, presented in Eq. 10, 
over the relativistic wave function for electrons near the 
nucleus 6] at zero energy. The use of the wave function 
for the electrons at zero energy is justified by noting that 
the main contribution to the matrix element of SH comes 
from distances close to the nucleus, r ~ a/Z. At a dis- 
tance of one Bohr radius, r ~ a, the potential is screened 
by the other electrons and the potential energy is given 
by Vc ~ — mc 2 o: 2 . This is of the same order of magni- 
tude as the binding energy of the electron, E ~ — mc 2 a 2 . 
However, close to the nucleus, at r ~ a/Z, screening 
is negligible and Vc ~ —Z 2 mc 2 a 2 . Inside this region 
\Vc\ ^S> \E\ and so the binding energy of the electron 
can be safely ignored. A detailed derivation of the Z 2 
enhancement of < SH > for r < a/Z is provided in the 
appendix. 

The non-relativistic limit of SE was taken. This gave: 



5E_ _ 2(Za) 2 tan 2 X 1 
~E ~ v~ 



3 + 1/2 



(4) 



where E = — Z 2 mc 2 a 2 / (2v 2 ) is the energy of the elec- 
tron, and Z a is the charge "seen" by the electron - it is 1 
for atoms, 2 for singly charged ions etc. The derivation 
of Eq. Q can be found in the appendix. 

It is interesting to note that this correction to the en- 
ergy of the electron has exactly the same form as the 
relativistic correction, A, to the energy of an external 
electron: 



A 

E 



(Za) 2 



1 



.7 + 1/2' 



(5) 



We can sum up Bekenstein's relativistic correction, Eq. 
(0J , and the usual relativistic correction to give 



~E 



(Za') 



IV2 



J + 1/2' 



(6) 



where a' — a(l+ tan 2 \)- 

The works 0, 0, Q used a method suggested in Q| for 
the analysis of absorption lines. A comparison between 



2 



different frequencies is used. In this method only the 
relativistic corrections, A' / E, are used to determine a 
variation since any variation in the energy in the non- 
relativistic limit is absorbed into the redshift parame- 
ter and also scales the same way for all elements. Since 
Eq. and 10 are directly proportional the effect of 
the modified form of the Dirac Hamiltonian is indistin- 
guishable from a small change in a 2 in Eq. JSJ. A 
measurable change in a would in fact be a change of 
a' = a(l + tan 2 x)- The astrophysical data can not dis- 
tinguish between a and tan 2 \ variation. 

Note that the proportionality of 5E and A has a 
simple explanation. The relativistic corrections to the 
Schrodinger Hamiltonian (e.g. the spin-orbit interaction) 
and the non-relativistic limit of SH have similar depen- 
dence on r and both are proportional to the electron 
density ip 2 at r ~ a/Z. The proportionality of SE and A 
also holds for high orbitals (small binding energy) in the 
pure Coulomb case (see next section). 

In this derivation we assumed that we could consider 
the unscreened Coulomb field, this is clearly not the case 
for a valence electron in a many electron atom. We justify 
this assumption by once again noting that the main con- 
tribution to SE is given by distances close to the nucleus, 
r ~ a/Z . At this distance the main screening comes from 
the Is electrons and we can use Slater's rules to estimate 
the screening corrections to Eq. (0J and JSJ: 



SE _ 2a 2 (Z-0.6) 2 tan 2 x 1 
~~E ~ v 



3 + 1/2 



A _ a 2 (Z-0.6) 2 
E ~ v 



1 



i + 1/2 



(7) 



(8) 



This does not affect the proportionality of the two terms 
and makes very little difference to the results in heavy 
atoms (~ 1/Z). The correction from the non-zero energy 
of the external electron is even smaller (~ 1/Z 2 ). Consid- 
eration of many-body correlation corrections has shown 
that this does not change the proportionality relationship 
either. The point is that the expressions for the corre- 
lation corrections obtained using the many-body pertur- 
bation theory (or the configuration interaction method) 
contain the single-particle matrix elements of SH and 
that of the relativistic corrections which are proportional 
to each other. This makes the final results proportional. 
Because of this proportionality it is not possible to derive 
values for a and tan 2 \ separately in mulitelectron atoms. 
However, a separation of these values can be achieved in 
hydrogen. 



III. CALCULATIONS INVOLVING HYDROGEN 
ATOM 

The case is somewhat simplified for the hydrogen atom 
and other hydrogen-like ions as there are no inter-electron 
interactions. There are also very accurate experimental 
measurements of transition frequencies in hydrogen. 



We confirm Bekenstein's result Q that applying the 
Hamiltonian Q one can derive: 



SE = - 



vac Z or 



^)tan 2 X . (9) 



n 3 Vi + 1/2 2n 
The relativistic correction to the electron energy is 



A = - 



2 4 

mc Z a 



1 



2n 3 \j + l/2 An 



(10) 



Note that for large n (zero energy), SE is again propor- 
tional to the relativistic correction, A. 

It is possible to obtain a limit on the tan 2 \ parameter 
by comparing the theoretical and experimental data for 
the 2p 3 / 2 — splitting. The experimental value 

/ 2p3/2 -2 Pl/2 (exp) = 10 969 045(15) kHz (11) 

is derived from two experimental results, 



/2s 1/2 ->2 P3/2 (exp) 
j2 Pl/2 ^2s 1/2 (exp) 



9 911 200(12) kHz (12) 
1057 845(9) kHz (13) 



presented in and Q respectively. It has to be com- 
pared with a theoretical value which we take from a com- 
pilation (see also review [Tof) 

/2 P3/2 -2p 1/2 (theory) = 10 969 041.2(15) kHz . (14) 

Noting that SE has not been accounted for in the Eq. 
(tT4"jl . but will be present in Eq. (fTTjl . and using Eq. © 
we can write: 

■ tan 2 x = /2 P3/2 ^2p 1/2 (exp)-/2p 3/2 ^2p 1/2 (theory). 



16h 



(15) 

The leading contribution to /2p 3/2 — >2pi /2 (theory) is given 
by mc 2 a 4 /32h, this allows us to write: 

/2p 3/2 -^2p 1/2 (exp) = / 2 p 3/2 ^ 2 p 1/2 (theory) (1 + 2tan 2 x) 

Using the values above we can obtain the limit tan 2 \ = 
2(7) x 10~ 7 . It is assumed here that /2p 3/2 -^2p 1/2 (theory) 
is expressed in terms of a which is extracted from the 
measurements of parameters which are not sensitive to 
tan 2 x (i- e - they depend on a rather than on a 1 ). In- 
deed, one of the values of a is derived via a complicated 
chain of relations with a eventually coming from the Ryd- 
berg constant which is quite weakly affected by SH (the 
relative value of the correction is of order of a 2 tan 2 x 
since the matrix element of SH vanishes in a leading 
non-relativistic approximation). The most accurate re- 
sult obtained this way is a' 1 = 137.036 000 3(10) 

A self consistent quantum electrodynamic theory with 
a dynamically varying a should meet some even stronger 
constraints due to a comparison of the value of the fine 
structure constant from the anomalous magnetic moment 
of the electron (a' 1 = 137.035 998 80(52) 12]) with the 
Rydberg constant value. Such a comparison will likely 



3 



lead to a limitation on tan 2 x at a level of a few parts in 
10~ 8 since 8a/ a — 11(8) x 10~ 9 , from comparison of the 
values for a^ 1 given above. 

Before any modification of QED due to a varying a 
is considered seriously another set of questions need to 
be answered. These questions should target its gauge 
invariance, renormalizability and Ward identities, which 
supports the same charge for electrons and protons. The 
current QED construction is quite fragile and it is not 
absolutely clear if it can be successfully extended. 



IV. CONCLUSIONS 

In conclusion, using the modified form of the Dirac 
Hamiltonian Eqs. Q~© does not affect the analy- 
sis used in 0, 0, 0. They measure the variation of 
a! = q(1 + tan 2 x). The present time limit tan 2 \ = 
(0.2 ± 0.7) x 10~ 6 is obtained from the measurement of 
the hydrogen 2p fine structure using value of a obtained 
from different experiments. Note that according to 
the value of a' was smaller in the past, the last measure- 
ment gave Aa'/a' = (-0.54 ±0.12) x 10~ 5 . If there is no 
other source of variation of a this would require a nega- 
tive value of tan 2 x (tan 2 x = (-0.52±0.14) x 10~ 5 ) since 
the present value of tan 2 \ is small. Actually, the choice 
of the integration constants in the Bekenstein paper pre- 
cludes considering epochs with a' < a 5]. However this 
should not be deemed a principle problem. 
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APPENDIX - CALCULATION OF SE 

The most convenient way to calculate SE is to calculate 
the matrix element of the operator SH = {1 — 0) tan 2 x ■ 
Vc for an external electron in a many-electron atom or 
ion using a relativistic wave function. In order to use the 
relativistic wave functions for electrons near the nucleus 
at zero energy it is necessary to demonstrate that the 
major contribution to the matrix element of SH comes 
from distances r < a/Z, where the screening of the nu- 
clear potential and the external electron energy can be 



neglected. Only the contribution at distances 



< 



a/Z 



have a Z 2 enhancement, the contribution at r ~ a does 
not have this enhancement since the atomic potential at 
this distance is screened, Vc ~ e 2 /r, and has no Z de- 
pendence. 

To demonstrate the Z 2 enhancement of the r < a/Z 
contribution let us consider the non-relativistic limit of 



the operator SH. The matrix 

1-0: 





-2 



has only lower components, it follows that the matrix 
element, 



tP + {I-0)V c iP = -2 X + VcX, 



(17) 



where 



is the Dirac spinor. In the non-relativistic limit 

<T ■ p 



X 



2mc 



<P 



and this gives 



SH 



1 



2m 2 c 2 



(CT-p)Vb(o--p) tan 2 x 



-Vc 



ih 



2m 2 c 



2m 2 c 



r(VVc-p) 



h 2 1 dV c 
2m 2 c 2 r dr 



a ■ 1 ) tan x- 



(18) 



This derivation is similar to the standard derivation of 
the spin-orbit interaction term in the non-relativistic ex- 
pansion of the Dirac Hamiltonian. Let us now compare 
the contributions of r < a/Z and r ~ a to the matrix 
element of SH. Consider, for example, the last spin-orbit 
term in Eq. I|18fl which is proportional to the usual spin- 
orbit interaction. The electron wave function at r ~ a/Z 
is given by ip 2 ~ Z/a 3 [13j. the spin-orbit operator is 
proportional to 



IdVc 



r dr 



Kze 2 



and the integration volume is proportional to r 3 . As a 
result we can write: 



2 C 2 



< SH >~ Z 2 



We 



(19) 



2m 2 c 2 a 3 

For r ~ a, the wave function is ip 2 ~ 1/a 3 , the spin-orbit 
operator is proportional to 



l dV c 
r dr 



1 



and the integration volume is still proportional to r 3 . 
Therefore 



2 C 2 



< SH >r 



2m 2 c 2 a 3 ' 



(20) 



this is Z 2 times smaller than at r < a/Z. The same con- 
clusion is also valid for the first two terms in Eq. JTHJl. 
This estimate demonstrates that the main contribution 
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to SH comes from small distances r < a/Z. This conclu- 
sion is similar to that for the relativistic corrections to 
atomic electron energy. 

We perform the actual calculation of the matrix ele- 
ment of SH using the relativistic Coulomb wave functions 
for zero-energy electrons near the nucleus. These can be 
expressed in terms of Bessel functions as |^ : 



fnjl(r) = ^[ja + k)J 2i {x) - |j 27 -i(») 
9njl(r) = —^-ZaJ 2y (x) 



(21) 



where x = (SZr/a) 1 / 2 , 7 = + 1/2) 2 - Z 2 a 2 , k 
(_l)J+i/2-iy + i/2) and 



1 \ 1/2 



|k| V Zav 3 



z a . 



SE can now be calculated using these wave functions and 
SH: 

SE = J 4i + 5HiPdV 

= -2tan 2 x / gnu v c9njir 2 dr 



r 

dx 



JO 

f°° dr 
= -2e 2 Z 3 a 2 tan 2 X c 2 njl / J 2 2 7 (x)- 
Jo 

/>oo 

= -4e 2 Z 3 « 2 tan 2 X c^ / J 2 ^)- 
Jo x 



We now use the relationships between Bessel functions 
and Gamma functions to write this as: 



■1 '/•").. ^ ,, p 2 ^(27) 
nj7 r(2 7 + i) 



SE = -2e 2 Z 3 a 2 tan 2 yc 2 



e 2 Z z a 2 tan 2 \c. 



rrjl 



When we substitute in for 



njl 



we obtain: 



SE = - 



mc 2 Z 2 Z 2 a tan 2 \ 



Finally, we take the non-relativistic limit by replacing 7 
with j + 1/2, 



SE = - 



mc z Z z Z z a a* t&n z x 
^ + 1/2) ■ 



(22) 



To obtain Eq. |@J wc simply divide Eq. by E = 

— Z 2 ma 2 /(2v 2 ). Note that the spin-orbit contribution 
to Eq. © can be obtained in an analogous manner. 

Finally, we should present a very simple derivation of 
Eqs. (0| and JSJ) based on the results obtained for the 
pure Coulomb case, see Eqs. (0 and (|10|l . For the high 
electron orbitals (n >> 1) the electron energy at r ~ 
a/Z may be neglected in comparison with the Coulomb 
potential and the Coulomb results Eqs. © and pO(l are 
proportional to the electron density at r ~ a/Z where 
■ijj 2 ~ 1/n 3 . For the external electron in heavy atoms the 
situation is similar. The external electron wave function 
in Eqs. (|21|l at r ~ a/Z is proportional to the Coulomb 
wave function for small energy (n >> 1). Therefore, 
to find the matrix elements for the external electron we 
should take the Coulomb results Eqs. © and l|10[l and 
multiply them by the ratio of the electron densities for 
the external electron in the many-electron atom and the 
Coulomb electron. This immediately gives Eqs. an d 

ta- 
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